't Hooft Determinant: Fluctuations and Multiple Vacua 
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The 't Hooft six quark flavor mixing interaction (TV/ = 3) is bosonized by the path integral 
method. The considered complete Lagrangian is constructed on the basis of the combined 't Hooft 
and U($) x !7(3) extended chiral four fermion Nambu - Jona-Lasinio interactions. The method of 
the steepest descents is used to derive the effective mesonic Lagrangian. Additionally to the known 
lowest order stationary phase (SP) result of Reinhardt and Alkofer we obtain the contribution from 
the small quantum fluctuations of bosonic configurations around their stationary phase trajectories. 
Fluctuations give rise to multivalued solutions of the gap equations, marked at instances by drastic 
changes in the quark condensates, a novel scenario for the vacuum state of hadrons at low energies. 



I. INTRODUCTION 



The global C/z,(3) x Ur(3) chiral symmetry of the QCD Lagrangian (for massless quarks) is broken by the £Za(1) 
Adler-Bell-Jackiw anomaly of the ST/ (3) singlet axial current q^^^q. Through the study of instantons [1,2], it became 
clear that effective 2Nf quark interactions, known as 't Hooft interactions, violate Ua(X), but are still invariant under 
\Q [ chiral SU(Nf) x SU(Nf). In the case of two flavors they are four-fermion interactions, and the resulting low-energy 
theory resembles the old Vaks Larkin Nambu Jona-Lasinio model [3] . In the case of three flavors they are six-fermion 
interactions which are responsible for the correct description of n and rl physics, and additionally lead to the OZI- 
violating effects [4,5], 



Alct = K(detqP R q + ActqP L q) (1) 



^ ' where the matrices Pr,l = (1 ± 7s)/2 are projectors and the determinant is over flavor indices. 
, The physical degrees of freedom of QCD at low-energies are mesons. The bosonization of the effective quark 
f^i ■ interaction (1) by the path integral approach has been considered in [6], where the lowest order stationary phase 
\ approximation (SPA) has been used to estimate the leading contribution from the 't Hooft determinant. In this 
approximation the functional integral is dominated by the stationary trajectories r s ^(x), determined by the extremum 
condition SS(r) = of the action S(r). The lowest order SPA corresponds to the case in which the integrals associated 
Qh, with S 2 S(r), for the path r s ^(x) are neglected and only S(r s ^) contributes to the generating functional. The subject 
of the present work is to include the contribution associated with 8 2 S(r). Our motivation is the following: let us 



consider the one-dimensional integrals of the form 

F(X) = I dtexp[iXf(t)} (2) 



What is sought by the method of SP is the dominant contribution to F as A — > oo. The dominant contribution to the 
integral comes from regions of t where /' vanishes. Supposing that /' vanishes at only one point to and neglecting 
contributions to the integral from regions far from to, one can obtain the result 



Hence F(X) goes to zero like 1/vA. The lowest order SPA would correspond to the result F(X) — const ■ exp[iA/(to)]j 
which has the incorrect asymptotic behavior F(X) — 0(1). If results related to finite-dimensional integrals, such as 
■FX A), mean anything with regard to corresponding functional integrals, one can conclude that the leading SPA should 
include the contribution from the integrals associated with the second functional derivative S 2 S(r s ^). 



*On leave from the Joint Institute for Nuclear Research, Laboratory of Nuclear Problems, 141980 Dubna, Moscow Region, 
Russia. 
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II. PATH INTEGRAL BOSONIZATION 



To be definite, let us consider the theory of the quark fields in four dimensional Minkowski space, with dynamics 
defined by the Lagrangian density 

£ = %JL + £def ( 4 ) 

The first term here is the extended version of the Nambu - Jona-Lasinio (NJL) Lagrangian £^ = £ + £ m t > 
consisting of the free field part 

C = qii^dft - rh)q, (5) 
and the U(3)l x U(3)r chiral symmetric four-quark interaction 

G 

Ant = j[(v x aq) 2 + [qilhKq?] ■ (6) 

We assume that quark fields have color and flavor indices running through the set i = 1,2,3; A a are the standard 
C/(3) Gell-Mann matrices with a = 0, 1, . . . , 8. The current quark mass, m, is a nondegenerate diagonal matrix with 
elements diag(m„, m<2, rh s ), it explicitly breaks the global chiral U(3)l x U(3)r symmetry of the Lagrangian. 
The second term in (4) is given by (1). Letting 

s a = -q\ a q, Pa = qhsKtQ, s = s a X a , p = p a \ a (7) 

yields 

£ det = I det ( s + *P) + dct ( s - ■ ( 8 ) 
The dynamics of the system is described by the vacuum transition amplitude in the form of the path integral 

Z = j VqVqexp (i J d^xC^j . (9) 

which can be written in the equivalent form [6] 



Z = y T><T a T>(j) a T>qT>qcxp J d^xCq^j j T>r\ a Vr2 a exp i^i J d i x£ r 



(10) 



where 



C q (q, q, a, cj>) = q{i^d^ - m - a + i-y^q, (11) 
C r (a,(t>,r 1 ,r 2 ) = 2G [(r la ) 2 + (r 2a ) 2 ] - 2(r la a a + r 2a 4>a) 

- 1 [det(n + ir 2 ) + det(n - ir 2 )] . (12) 
8 

The Fermi fields enter the action bilinearly, we can always integrate over them. At this stage one should also shift the 
scalar fields a a — > cr a + A a by demanding that the vacuum expectation values of the shifted fields vanish < 0|cr Q |0 >= 0, 
yielding gap equations to fix parameters A a = m a — m a , where m a denotes the constituent quark masses . To evaluate 
path integrals over 7*1,2 one nas to use the method of stationary phase, or, after the formal analytic continuation in 
the time coordinate X4 = ix , the method of steepest descents. The analytical continuation of the Euclidean version 
of the path integral under consideration is (see [8] for details), 

Vr la Vr 2a exp I / d 4 x£ r ((j, <j>, 7*1 , r 2 ) J . (13) 



'The shift by the current quark mass is needed to hit the correct vacuum state, see e.g. [7]. 
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Near the saddle point r^, 

C r « C r (r st ) + \ J2^C (r at )r f3 (14) 



2 

Q,/3 



where is a solution of the equations £J,(ri, r2) = determining a flat spot of the surface C r {r\,r2), 

/ 2Grf - (a + A) a - f A a6c (r*rJ - r£r§) =0 

\ 2Gr? - 4> a + ^A abc r\rl = 0. [10) 

with totally symmetric constants, A a b c , closely related with the U(3) constants d a bc- We use in (14) symbols f a for 
the differences (r a — r gt)- To deal with the multitude of integrals in (13) we define a column vector f with eighteen 
components r a = (f^ff) and with the matrix C'^(r s ^) being equal to 

£%{r s t) = *GQ a0 , Q a(i = ( 5ab ~ f fcjfff ) . (16) 

Eq.(13) can now be concisely written as 

J(a,</>) = exp ^ d^Crir^ J + *°° Vf a exp (^2G J d^x^Qir^ [1 + 0(h)] . (17) 

The solutions of Eqs.(15) are the following even and odd parity combinations and r° gt expressed in the form 
of increasing powers in a a , 4> a 

rlst = h « + h ab <*> + h al^c + h^Mc + ■■■ (18) 

^St= /l i6^ + ^ 3 6 ) c^c + ... (19) 

Putting these expansions in Eqs.(15) one obtains a series of selfconsistent equations to determine the constants h a , 
hail > h a 2 b \ etc. [8] As a result we get 

£ r (r st ) = -2h a a a - h^a a a b - <j> a <j>b + ©(field 3 ). (20) 

We now turn to the evaluation of the path integral in Eq.(17). In order to define the measure Vf a more accurately 
we expand f a in a Fourier series 

oo 

r<x{x) = Cn,atPn(x), (21) 
n=l 

assuming that suitable boundary conditions are imposed. The set of the real functions {f n (x)} form an orthonormal 
and complete sequence, therefore 

fvf a exp Ug f d^Q^f] = ° (22) 

J \ J / v /det(2GA^) 

2G\f m = ( ^ 1 fc(2 °-i J (s <r pS nm + J d*xM*)F°p{*)v m (x)) (23) 



with 

(L, 1 



and 

F„-^A H (- h -^U- K) r!® 1 ** ) (24) 

-tV/3 — . Sieba (2) , (2)/ „ , \ ' 

Only the matrix F CT/ g depends here on fields a, 4>. By absorbing in C the irrelevant field independent part of 2GX^ nl 
and expanding the logarithm in the representation det(l + F) = exptr ln(l + F), one can obtain finally for the complete 
action 

Sr = J d^x |c(r st ) + ^ £ [F^(r st )] j (25) 

proposing that the undetermined dimcnsionlcss constant a will be fixed by confronting the model with experiment 
afterwards [9] 
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III. THE GROUND STATE 



Let us study the ground state of the model under consideration, then properties of the excitations will follow 
naturally. In Eq.(12) the field coefficients hi obey 

2Ghi - Aj = -Ujkhjhk (26) 

o 

with the totally symmetric coefficients t^ k equal to zero except for t uds = 1 and with order parameters A, ^ (i = 
u, d, s). 

A tadpole graphs calculation gives for the gap equations the following result 

2fc + g (*£ - A*,^ - (27) 

where the left hand side is the contribution from (25) and the right hand side is the contribution of the quark loop 
from (10), [8]. 

The second term on the left hand side of Eq.(27) is the correction resulting from the Gaussian integrals of the 
steepest descent method, comprising the effects of small fluctuations around the stationary path. If one puts for a 
moment a = in Eq.(27), and combines the result with Eqs.(26), one finds gap equations which are very similar to 
the ones obtained in [4]. For this case one obtains values of (m u , m s ) which arc uniquely related to values of (G, n). 

The general case, which we have when a > in Eq.(27), yields in turn a region for m u ,m s , where three values of 
couplings (G, n) are possible. 

Conversely, one can study the solutions: m u = m u (G,K), m s = m s (G,K) at fixed values of input parameters: 
A, m u = rhd, rh s . Again we find several extremal solutions. The minima will be identified in a future analysis of the 
effective potential. We refer to [8] for graphical displays. 

Multivalued solutions of the gap equations have been obtained in a different context in [10]. 



IV. CONCLUDING REMARKS 



The purpose of this work has been twofold. Firstly we have developed the technique which is necessary to go beyound 
the lowest order SPA in the problem of the path integral bosonization of the 't Hooft six quark interaction. This 
technique is rather general and can be readily used in other applications. Second, we have explored with considerable 
detail the implications of taking the quantum fluctuations in account in the description of the hadronic vacuum. We 
encountered several classes of solutions. These multiple vacua may have very interesting physical consequences and 
applications. 
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